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Abstract 

In this series of eight papers we present the applications of 
methods from wavelet analysis to polynomial approxima- 
tions for a number of accelerator physics problems. In this 
part we consider a model for spin-orbital motion: orbital 
dynamics and Thomas-BMT equations for classical spin 
vector. We represent the solution of this dynamical sys- 
tem in framework of biorthogonal wavelets via variational 
approach. We consider a different variational approach, 
which is applied to each scale. 

1 INTRODUCTION 

This is the fourth part of our eight presentations in which 
we consider applications of methods from wavelet analy- 
sis to nonlinear accelerator physics problems. This is a 
continuation of our results from [l]-[8], which is based on 
our approach to investigation of nonlinear problems - gen- 
eral, with additional structures (Hamiltonian, symplectic or 
quasicomplex), chaotic, quasiclassical, quantum, which are 
considered in the framework of local (nonlinear) Fourier 
analysis, or wavelet analysis. Wavelet analysis is a rela- 
tively novel set of mathematical methods, which gives us 
a possibility to work with well-localized bases in func- 
tional spaces and with the general type of operators (dif- 
ferential, integral, pseudodifferential) in such bases. In this 
part we consider spin orbital motion. In section 3 we con- 
sider generalization of our approach from part 1 to varia- 
tional formulation in the biorthogonal bases of compactly 
supported wavelets. In section 4 we consider the different 
variational multiresolution approach which gives us possi- 
bility for computations in each scale separately. 

2 SPIN-ORBITAL MOTION 

Let us consider the system of equations for orbital motion 
and Thomas-BMT equation for classical spin vector [9]: 

dq/dt = dH or i,/dp, dp/dt = —dH or b/dq, ds/dt = 
w x s, where 
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q = 92, 93), P = (pi,P2,P3) are canonical position 
and momentum, s = (s ly s 2 , S3) is the classical spin vec- 
tor of length Ti/2, it — 7T2, W3) is kinetic momen- 
tum vector. We may introduce in 9-dimensional phase 
space z = (q,p,s) the Poisson brackets {f(z),g(z)} = 
fq9p ~ fp9q + [fs x 9s] ■ s and the Hamiltonian equations 
are dz/dt = {z, H} with Hamiltonian 



H = H ar b(q,p,t) + w(q,p,t) ■ s. 
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More explicitly we have 
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dp 
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We will consider this dynamical system also in another pa- 
per via invariant approach, based on consideration of Lie- 
Poison structures on semidirect products. But from the 
point of view which we used in this paper we may con- 
sider the similar approximations as in the preceding parts 
and then we also arrive to some type of polynomial dynam- 
ics. 

3 VARIATIONAL APPROACH IN 
BIORTHOGONAL WAVELET BASES 

Because integrand of variational functionals is represented 
by bilinear form (scalar product) it seems more reasonable 
to consider wavelet constructions [10] which take into ac- 
count all advantages of this structure. The action functional 
for loops in the phase space is [11] 



H or b 
w = 



cy/ ir 2 + 



e<I> 



TO0C7 



-(1 + 7 G)£ 



(1) 



* e-mail: zeitlin@math.ipme.ru 

t http://www.ipme.ru/zeitlin.html; http://www.ipme.nw.ru/zeitlin.html 



F( 7 ) 



pdq 



H(t,-y(t))dt 



(4) 



The critical points of F are those loops 7, which solve 
the Hamiltonian equations associated with the Hamiltonian 



H and hence are periodic orbits. By the way, all critical 
points of F are the saddle points of infinite Morse index, 
but surprisingly this approach is very effective. This will 
be demonstrated using several variational techniques start- 
ing from minimax due to Rabinowitz and ending with Floer 
homology. So, (M, ui) is symplectic manifolds, H : M — > 
R, H is Hamiltonian, Xh is unique Hamiltonian vec- 
tor field defined by u>(Xh(x), v) — —dH(x)(v), v G 
T X M, x G M, where u> is the symplectic structure. 
A T-periodic solution x(t) of the Hamiltonian equations 
x = Xh(x) on M is a solution, satisfying the boundary 
conditions x(T) = x(0),T > 0. Let us consider the loop 
space n = C 0O (5 1 ,i? 2 ™), where S 1 = R/Z, of smooth 
loops in R 2n . Let us define a function $ : fl — > R by 
setting 

= / ^ < - Jx, x > dt- [ H(x(t))dt, x ett 
Jo 2 Jo 

(5) 

The critical points of $ are the periodic solutions of x = 
X H (x). Computing the derivative at x e O in the direction 
of y G f2, we find 



$'(*)(!/) - ^(x + ey)\ £=0 
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< - Jx - SjH(x), y > dt 



Consequently, = for all y e O iff the loop x 

satisfies the equation 



Jx{t) - S/H{x(t)) = 0, 



(7) 



i.e. x{t) is a solution of the Hamiltonian equations, which 
also satisfies x(0) = x(l), i.e. periodic of period 1. Pe- 
riodic loops may be represented by their Fourier series: 
x(t) = e fc27rJ *Xfe, Xk G R 2k , where J is quasicomplex 
structure. We give relations between quasicomplex struc- 
ture and wavelets in our other paper. But now we need to 
take into account underlying bilinear structure via wavelets. 
We started with two hierarchical sequences of approxima- 
tions spaces [10]: 

...V- 2 CV- 1 CV CV 1 CV 2 ..., (8) 
..T_ 2 cy_ickVici., 

and as usually, Wo is complement to Vb in V\, but now 
not necessarily orthogonal complement. New orthogonal- 
ity conditions have now the following form: 

Wo -L Vb, Wo -L V , Vj 1 Wj, Vj 1 Wj (9) 



where ipjk{x) — 2 : >/ 2 ip{2^ x — k). Functions if(x), (p{x—k) 
form dual pair: < ip(x — k),<p(x — £) >= Ski, < <p(x — 
k),ip(x — t) >= 0. Functions ip, (p generate a multiresolu- 
tion analysis. (p(x — k), ip(x — k) are synthesis functions, 
(p(x — £), tp(x — £) are analysis functions. Synthesis func- 
tions are biorthogonal to analysis functions. Scaling spaces 
are orthogonal to dual wavelet spaces. Two multiresolu- 
tions are intertwining Vj+Wj = Vj+i, Vj+Wj = Vj+i. 
These are direct sums but not orthogonal sums. 

So, our representation for solution has now the form 



/(*) = X^' fe ^' fc W' 



(11) 



where synthesis wavelets are used to synthesize the func- 
tion. But bjk come from inner products with analysis 
wavelets. Biorthogonality yields 
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f{t)i> tm (t)dt. 
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So, now we can introduce this more complicated construc- 
tion into our variational approach. We have modification 
only on the level of computing coefficients of reduced non- 
linear algebraical system. This new construction is more 
flexible. Biorthogonal point of view is more stable un- 
der the action of large class of operators while orthogonal 
(one scale for multiresolution) is fragile, all computations 
are much more simpler and we accelerate the rate of con- 
vergence. In all types of Hamiltonian calculation, which 
are based on some bilinear structures (symplectic or Pois- 
sonian structures, bilinear form of integrand in variational 
integral) this framework leads to greater success. In par- 
ticular cases we may use very useful wavelet packets from 
Fig. 1. 




Figure 1: Wavelet packets. 



translates of tp span Wo, translates of tp span Wo- 
Biorthogonality conditions are 

/oo 
ipjk(x)ip j 'k'(x)dx = Skk'Sjj', 
-OO 

(10) 



4 EVALUATION OF NONLINEARITIES 
SCALE BY SCALE.NON-REGULAR 
APPROXIMATION. 

We use wavelet function tp(x), which has k vanishing mo- 
ments J x k ip(x)dx — 0, or equivalently x k = ^2c(<p((x) 



for each k, < k < K. Let Pj be orthogonal projector on 
space Vj . By tree algorithm we have for any u G L 2 (R) 
and f £ Z, that the wavelet coefficients of P&(u), i.e. the 
set {< u, ipj^k >,j < i — 1, k € Z} can be compute using 
hierarchic algorithms from the set of scaling coefficients in 
Vi, i.e. the set {< u,ifi£^ >,k G Z} [12]. Because for 
scaling function y> we have in general only J ip(x)dx = 1, 
therefore we have for any function u G i 2 (R): 

2 o/2 < u,ip jyk > -u(x) |= 
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If the integer n(tp) is the largest one such that 



x a <p(x)dx = 



for 



1 < a < n 



(13) 
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then if u £ C^ n+1 ' ) with u( n+1 ) bounded we have for j — > 
oo uniformly in k: 

| 2 J/2 < u, <ft, fc > -u(A;2- J ') |= 0(2^' ( " +1) ). (15) 

Such scaling functions with zero moments are very use- 
ful for us from the point of view of time-frequency local- 
ization, because we have for Fourier component $(w) of 
them, that exists some C(ip) £ R, such that for uj — > 
$(w) = 1 + C(ip) | cj | 2r + 2 (remember, that we con- 
sider r-regular multiresolution analysis). Using such type 
of scaling functions lead to superconvergence properties 
for general Galerkin approximation [12]. Now we need 
some estimates in each scale for non-linear terms of type 
u i — > f(u) = f o u, where f is C°° (in previous and future 
parts we consider only truncated Taylor series action). Let 
us consider non regular space of approximation V of the 
form 

V = V q ® W J> < 16) 

q<j<p-l 

with Wj C Wj . We need efficient and precise estimate of 
/ o u on V. Let us set for q G Z and u G L 2 (R) 

Y[ /,(«) = 2-^ 2 f(2 q/2 < u, Vq , k >) ■ ^ k (17) 
feez 

We have the following important for us estimation (uni- 
formly in q) for u, f(u) G [12]: 



\\P 9 (/(«)) -n/^)IU 2 = °( 2 ~ ( " +1)< " 8 » 
For non regular spaces (|l6[) we set 
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Then we have the following estimate: 

||P^ (/(«)) -Il/yWIU 2 =0(2-("+^) (20) 

uniformly in q and V (|l6|). This estimate depends on q, not 
p, i.e. on the scale of the coarse grid, not on the finest grid 
used in definition of V. We have for total error 

+ \\Py(f(u)-Y[f^(u))\\ L2 (21) 



and since the projection error in V: \\f(u) — P v (f(u)) \\^2 
is much smaller than the projection error in V q we have 
the improvement (20) of (|l8|). In concrete calculations and 
estimates it is very useful to consider approximations in the 
particular case of c-structured space: 



V = V Q + spanii/j-j.k, 



(22) 
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